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Abstract 



In this paper, we will define a quantum operator that performs the 
inversion about the mean only on a subspace of the system {Partial Dif- 
fusion Operator). This operator is used in a quantum search algorithm 
j^— ^ ' that runs in 0{\/N/M) for searching an unstructured list of size A'^ with 

^b , M matches such that 1 < M < A'^. We will show that the performance of 

f~^ ' the algorithm is more reliable than known fixed operators quantum search 

algorithms especially for multiple matches where we can get a solution af- 
ter a single iteration with probability over 90% if the number of matches 
is approximately more than one-third of the search space. We will show 
_i ^ that the algorithm will be able to handle the case where the number of 

rS . matches M is unknown in advance such that 1 < M < Af in 0{\/ N/M). 

r—{ ' A performance comparison with Grover's algorithm will be provided. 

^ 1 Introduction 

}|J^ • In 1996, Lov Grover [12] presented an algorithm for searching an unstructured 

5t ! list of N items for a single match vifith quadratic speed-up over classical algo- 

rithms. Grover's original algorithm exploits quantum parallelism by preparing 
a uniform superposition that represents all the items in the list then iterates 
both an oracle that marks the desired item by applying a phase shift of -1 on 
that item (e-^^ , with 6i = tt) and nothing on the other items (e-^^ , with 02 = 0) 
and an operator that performs inversion about the mean (diffusion operator) to 
amplify the amplitude of the match, the processes of this operator includes the 
operation (2 |0) (0| — /) which applies a phase shift of -1 on the states within 
the superposition {e-'^^ , with 0i ~ n) except the state |0) " where it applies 
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nothing (e-*^^ , with (/)2 = 0) ^Sl ■ To maintain consistency with hteratures, this 
operation can also be written as (/ — 2 |0) (0|) which apphes a phase shift if -1 
on the state |0) (e-'^^, with 02 = tt) and nothing on the other states of the 
superposition {e-'^^ , with 0i = 0) together with a global phase shift of -1 |15j . 

It was shown that the required number of iterations is approximately 7r/4 yiV 
which is proved to be optimal to get the highest probability with the minimum 
number of iterations [201, such that there is only one match in the search space. 

In 1 181 1151 111! 1161 [2j, Grover's algorithm is generalised by showing that 
the uniform superposition can be replaced by almost any arbitrary superpo- 
sition and the phase shifts applied by the oracle and the diffusion operator 
(e-^i , e-^^ , e-"^^ and e-"^^ ) can be generalised to deal with the arbitrary superpo- 
sition and/or to increase the probability of success even with a factor increase 
in the number of iterations to still run in 0{yjN/M). These give a larger class 
of algorithms for amplitude amplification using variable operators from which 
Grover's algorithm was shown to be a special case. 

In another direction, work has been done trying to generalise Grover's algo- 
rithm with a uniform superposition for the case where there are known number 
of multiple matches in the search space ^ El , where it was shown that 
the required number of iterations is approximately tt /A^JN/M for small M/N. 
The required number of iterations will increase for M > A^/2, i.e. the problem 
will be harder where it might be excepted to be easier ^Hl- Another work has 
been done for known number of multiple matches with arbitrary superposition 
and phase shifts ^1 121 El ^^ E] where the same problem for multiple matches 
occurs. In El Eli a hybrid algorithm was presented to deal with this prob- 
lem by applying Grover's fixed operators algorithm for -k j^^ N jM times then 
apply one more step using different oracle and diffusion operator by replacing 
the standard phase shifts 6'i = tt, 0i = tt with accurately calculated phase shifts 
Q^ and <\)-^ according to the knowledge of the number of matches M to get the 
solution with probability close to certainty. Using this algorithm will increase 
the hardware cost since we have to build one more oracle and one more diffu- 
sion operator for each particular M. For the sake of practicality, the operators 
should be fixed for any given M and are able to handle the problem with high 
probability whether or not M is known in advance. 

In the case of multiple matches, where the number of matches is unknown, 
an algorithm for estimating the number of matches (known as quantum counting 
algorithm) was presented JHIEI- In Bl? another algorithm was presented to find 
a match even if the number of matches is unknown which will be able to work 
if M lies within the range 1 < M < 3iV/4. 

In this paper, we will propose a fixed operator quantum search algorithm 
that is able to handle the whole range 1 < M < A^ more reliably whether 
or not the number of matches in known in advance. The plan of the paper 
is as follows: Section 2 introduces the general definition of the unstructured 
search problem. Section 3 defines the partial diffusion operator [2]. Section 4 
introduces the algorithm and an analysis on its behaviour. Section 5 shows a 
comparison with Grover's algorithm. Section 6 introduces the algorithm shown 
in 0] for unknown number of matches by replacing Grover's algorithm with the 



algorithm proposed here. The paper wiU end up with a general conclusion in 
Section 7. 



2 Unstructured Search Problem 

Consider an unstructured list L oi N items. For simplicity and without loss of 
generality we will assume that N — 2"" for some positive integer n. Suppose 
the items in the list are labelled with the integers {0, 1, ...,N — 1}, and consider 
a function (oracle) / which maps an item z g L to either or 1 according to 
some properties this item should satisfy, i.e. f : L ^ {Ojl}- The problem is 
to find any i £ L such that f{i) = 1 assuming that such i exists in the list. In 
conventional computers, solving this problem needs O (N/M) calls to the oracle 
(query), where M is the number of items that satisfy the oracle. 

3 Partial Diffusion 

In this section, we will define the Partial Diffusion Operator Dp which performs 
the inversion about the mean only on a subspace of the system. The diagonal 
representation of the partial diffusion operator Dp when applied on n + 1 qubits 
system can take this form: 

Dp = (iJ®" ® h) (2 |0) (0| - /„+i) (if®" ® h) , (1) 

where the vector |0) used in Eqn. ^is of length 2N = 2"+^, Ik is the identity 

matrix of size 2*"' x 2*^ and H is the Hadamard gate i H — -y= j . 

To understand the effect of this operator, consider a general state \ip) of n+ 1 
qubits register: 

2JV-1 

H')=Y.5k\k). (2) 

fc=0 

For our purposes and without loss of generality, the general system \^) can 
be re-written as, 

N-l N-1 

iV') = E"^(i^')®|o)) + E/?^(ij')®|i))' (3) 

3=0 ]=0 

where {aj = 5^ : k even} and {/?-,- = 5^ : k odd}. The effect of applying Dp on 
\xl}) produces. 
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Figure 1: Quantum circuit representing the partial diffusion operator Dp over 
n + 1 qubits. 
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where (a) = ^7=0 '^i 1^ ^^ ^^*^ mean of the amplitudes of the subspace entan- 
gled with the state |0) of the extra qubit workspace, i.e. aj (|j) ® |0)). That 
is, applying the operator Dp will perform the inversion about the mean only 
on the subspace aj (|j) (81 |0)) and will only change the sign of the amplitudes 
for the rest of the system (3j {\j) ®\1)). A circuit implementation for Dp using 
elementary gates Q is shown in Fig. ^ 



4 The Algorithm 

In this section we will propose the algorithm assuming that the number of 
matches is known in advance. For a list of size N = 2^ , prepare a quantum 
register of size n + 1 qubits all in state |0) and apply the steps of the algorithm 
as follows (its quantum circuit is as shown in Fig. [SJ: 

1- Apply Hadamard gate on each of the first n qubits. 

2- Iterate the following steps q times: 



Apply the oracle Uj. 
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Figure 2: Quantum circuit for the proposed algorithm. 



ii- Apply the partial diffusion operator Dp. 
3- Measure the first n qubits. 

4.1 Analysis of Performance 

For the sake of clarity and to understand the behaviour of the algorithm, 

we will trace the algorithm during the first few iterations. The mechanism of 

amplifying the amplitudes can be understood as shown in Fig. |3| Now consider 

the algorithm if iterated once. Its behaviour can be understood as follows: 

1- Register Preparation. Prepare a quantum register of n + I qubits all in 

state |0), where the extra qubit is used as a workspace for evaluating the 



oracle Uf , the state of the system 



W, 



(i; 



can be written as follows, where 



the subscript number refers to the step within the iteration and (1) in the 
superscript refers to the iteration number: 



Wr 



f'>)^\Of''<g>\0). 



(5) 



2- Register Initialisation. Apply Hadamard gate on each of the first n qubits 
in parallel, so they contain the 2" states representing the list, where i is 
the integer representation of items in the list: 



N-l 



wl' 



H 



(g)n , 



w, 



(1) 



E I*) ® io) 



(6) 



3- Applying the Oracle. Apply the oracle Uf that maps the items in the list 
to either or 1 simultaneously and stores the result in the extra workspace 
qubit: 
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Figure 3: Mechanism of amplitude amplification for the proposed algorithm 
with N ^4 and A/ = 1. 
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4- Partial Diffusion. Apply the partial diffusion operator defined above. Let 
M be the number of matches, which make the oracle Uf evaluate to 1 
(solutions) such that 1 < M < N. Assume that ^^ ' denotes a sum over 
i which are desired matches, and ^^ " denotes a sum over i which are 



undesired items in the list. So, the system 
be written as follows: 



W2 ) shown in Eqn. [3 can 



JV-l 



N-l 



W. 



(1) 



^E"(N)®|o)) + ^E'(K)®|i)) 



(8) 



.(i; 



Applying Dp on W2 ) will result in a new system described as follows: 



N-1 N-1 N-1 



where the mean used in the definition of partial diffusion operator is, 



and ai, 61 and ci used in Eqn. El arc calculated as follows: 

«i = 2(ai)--^, 6i = 2(ai), ci = ^. (11) 

Such that, 



{N ~ M)al + Mbl + Mci^l. (12) 

Notice that, the states with amplitude 61 had amplitude zero before ap- 
plying Dp as shown in the first iteration in Fig. |S1 

5- Measurement. If wc measure the first n qubits after the first iteration 
{q = 1), then the probabilities of the system will be as follows: 

i- Probability Ps to find a match out of the M possible matches 
is given by taking into account that a solution \i) occurs twice as 
{\i) (g) |0)) with amplitude bi and {\i) |1)) with amplitude ci as 
shown in Eqn. ^ 



P^'^^M{bl + 4) 

= 5(#)-8(f)' + 4(^)^ 

ii- Probability Pns to find undesired result out of the states is given by: 

Pi'J = {N-M)al (14) 

Notice that, using Eqn. EJ 

PW+PW = 1. (15) 



Consider the system after first iteration shown in Eqn. before applying 
the measurement, the second iteration wiU modify the system as follows: 

Applying the oracle Uf will swap the amplitudes of the states which represent 
the matches, i.e. states with amplitudes bi will be with amplitudes ci and states 
with amplitudes ci will be with amplitudes bi so the system can be described 
as. 
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w^P) = «i E " (I*) ® io)) + ^1 E ' (I*) ® io)) + ^1 E ' (I*) ® ii))- (16) 
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Applying the operator Dp will change the system as follows. 



Wr 



(2) 
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«2E"(i*) 



N-l 
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^'2E'(i*)®|o))+^2E'(i*)®|i))' (17) 



1=0 i=0 4=0 

where the mean used in the definition of partial diffusion operator is, 

(a2) = ^((iV-M)ai + Mci), 
and a2, &2 a-nd C2 used in Eqn. 1171 arc calculated as follows: 



(18) 



a2 = 2(a2)-ai, 62 = 2(a2)-ci, C2 = -6i, 
and the probabilities of the system are, 



(19) 



P, 



(2) _ 



M {bl + cl) 
M{bl + bl), 



and, 



P^n'J ^(N- M) {al) 

= {N -M){b2+C2) 



(20) 



(21) 



= {N-M){b2-b,f. 
In the same fashion, the third iteration will give the following system. 



w['^)^Uf 



Wr 



(2) 



N-l N-l N-l (22) 

^2 E " (N) ® |o)) + C2 E ' (N) ® |o)) + 62 E ' (N) ® |i))- 

i=0 4=0 4=0 



W. 



(3) 



i:". 



M^f)) 



JV-l N-l N-l (23) 

as E "{\^®m+h E '(K)®|0))+C3 E '(N)®|1)), 

4=0 4=0 4=0 



where the mean used in Dr, is, 



(a3)-^((iV-Af)a2 + Mc2), (24) 

and as, 63 and C3 used in Eqn. 1231 arc calculated as follows: 

a3 = 2(a3)-a2, 63 = 2(a3)-C2, 03^-62, (25) 

and the probabilities of the system are, 



and, 



= M{bl + bl), 
Pi^) = (iv - M) (ai) 



(26) 



^{N- M) (63 + C3)' (27) 

{N - M) (b-s - 62)' . 



In general, the system after q > 2 iterations can be described using the 
following recurrence relations, 

JV-l N-l Af-1 

W^'^) = «. E "(I*) ® 10))+^'^ E '(I*) ® |0)) +^<^ E '(I*) ® 11))' (28) 

i=0 i=0 1=0 

where the mean to be used in the definition of the partial diffusion operator is 
as follows: Let y = 1 — M/N and s = \/^/N , then. 

{uq) =yaq-i + {l-y)cq-i, (29) 

and Cg, bq and c^ used in Eqn. 1281 are calculated as follows: 

ao = s, ai = s (2j/ - 1) , a, = 2 (a,) - a,_i, (30) 

5o = s, 6i = 2sy, 6, = 2 (a,) - Cq_i, (31) 

Co = 0, ci = -s, Cg = -6g-i, (32) 

and the probabilities of the system are. 



Pi'^)=A/(fe2 + c^) 



(33) 



and 

pi^) = (iV-Af)(a2) 

= (iV-M)(6, + c,)' (34) 

= (iV-Af)(&5-6,_i)'. 

Solving the above recurrence relations for a^, hq and c^ shown Eqn. I3UI Eqn. 
1311 and Eqn. |211 respectively, the closed forms are as follows: 

aq = s {Uq{y) ~ Uq^iiy)) , bq = sUq{y), Cq = -sUq-i{y), (35) 

where y = cos (6*) = 1 - M/N, < 9 < tt/2 and Uq{y) is the Chebyshev 
polynomial of the second kind JUj ^which is defined as follows. 



sin ((g + l ) 9) 
sin (9) 
The probabilities of the system. 



u<> (y) - 'IJ^:' ■ (36) 



P^'^^ = (1 - cos {9)) {U^{y) + U^_,{y)) , (37) 

and, 

P^f ^ cos (9) {Uq{y)-Uq^,{y))\ (38) 



^To clear ambiguity, Uf represents the oracle function and Uq{y) is the Chebyshev poly- 
nomial of the second kind. 
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Such that, 



P^'^^ + Pif = M (62 + cl) + {N-M) al 

= N{hl+ Cl) +2{N-M) Cgbg 

= idm (sin' {{q +1)0)+ sin' {qd) - 2 cos {9) sm{{q + 1) 9) sin {q9)) 
= ihlV) ('=°'' ('^^) '^^' (^) " '^^' ('?^) ^°'' (^) + '^"' ('^^)) 

_ sin'' (9) _ -, 

~ sin^(e) ~ -•-■ 

(39) 
Now, we have to calculate how many iterations, 5, are required to find any 
match with probability close to certainty for different cases of 1 < M < TV. To 
find a match with high probability on any measurement, then Ps must be as 
close as possible to one. To calculate the number of iterations, q, required to 
satisfy this condition, we need the following theorem. 

Theorem 4.1. Consider the following relation, 

Pf^'' = (1 - cos {9)) [{Ug {y)f + ([/^_i {y)f) = 1, (40) 

where Uq (jj) is the Chebyshev polynomial of the second kind, y — cos (9) and 
<9 < tt/2, then, 

- ^^^ a ^ 
q = or 9 ~ —. 

^29 2 

Proof. From the definition of Uq (y) shown in Eqn. 1361 then Eqn. I4UI can take 
this form. 



(1 _ cos (9) ) ( °"\^^;;^'^^ + ^^;v^ ) = 1, 



'sin2((5+l)6i) sin2(g6') 
sin^ (9) ^ sin2(6') 
or, 

sin^ {{q + 1)9)+ sin^ {q9) = 1 + cos (9) . 

Using simple trigonometric identities, the above relation may take the form, 

cos {2q9 + 29) + cos {2q9) + 2 cos (9) = 0. 

Using the addition formulas for cosine we get, 

2 cos {2q9) cos^ {9) ~ 2 cos {9) sin (2q9) sin (9) + 2 cos {9) = 0, 

2 cos {9) (cos {2q9) cos (9) - sin [2^9) sin [9) + 1) = 0, 

cos {9) (cos {2q9 + 6*) + 1) = 0. 

11 



From the last equation we get, 

cos (6*) = or cos(2g6' + 6') = -1, 

which gives the required conditions, 9 = 7r/2 or ^ = (tt — 0)/26. 

D 

The number of iterations must be integer, let q = \tt/26\ where |(Z — g| < 1/2. 
And since, cos (6*) = 1 - M/N, we have 9 > sin (6*) = V2NM - N'P/N, then, 

(41) 



TT 

29. 




(q) 



where [ J is the floor operation. To determine the lower bound of Pi using g, 
let Ps to take the form, 

P^'''='-^^{sinmq+l)9) + sin^q9)) 

= T+^{l~cosi9)cos{2qe + 9)). ^ ' 



We have, 



then. 



\q-q\< ^ 



2' 



\{2q + I) 9 - {2q + 1) 9\ <9, 
and from the definition of q, 

i2q+l)9^7r, 

then, 

cos((2q + l)6'-7r) < cos (6*) , 

or, 

-cos ((25+1)6*) <cos(6'). 

Using this in Eqn021 we get the following lower bound, 

p(,) > 1 + C0SH9) l + (l-#)^ .^3 

"^^ - l + cos(0) l + (l-#)- ^ ^ 

The minimum of the lower bound is 2V2 - 2 (« 0.83) when M/N = 2 - V2 
(« 0.5857). Notice that, when M/N w 0.5857, the probability of success is 
98.78% after a single iteration using Eqn. E| This minimum of the lower 
bound can be neglected with respect to the real behaviour of the algorithm. 

To demonstrate the real behaviour of the algorithm, we may plot the proba- 
bility of success Ps using the required number of iterations q for any given M, 
Fig. 21 shows this behaviour as a function of < M/N < 1. We can see from 
the plot that the minimum probability that the algorithm may reach is approxi- 
mately 87.88% when M/N « 0.2928. For 0.2928 < M/N < 1, g = 1 where only 
a single iteration is required to handle this range. For M/N < 0.2928, q > 1 
where the algorithm will behave more reliable than Grover's algorithm as we 
will see. 
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Figure 4: Probability of success of the algorithm using the required number of 
iterations q. 



5 Comparison with Grover's Algorithm 

First we will summarise the probabilities of success and the required number of 
iterations for Grover's algorithm and the proposed algorithm before giving the 
comparison. The probability of success of Grover's algorithm as shown in ^ is 
as follows: 



Pi«'^)=sin2((2gG + l)eG)>l 



M 

iV' 



(44) 



where sin^ (Og) = M/N, < 6'g < 7r/2 and the required number of iterations 
qa is, 



9G 



40G 



TT N 

- 4\Im' 



For the proposed algorithm, the probability of success is as follows, 



Pi^) = (l-cos(0)) 



sm^{{q+l)e) sm^iqd) 



+ 



sin^ (0) sin^ (9) 

where cos (9) = I — M/N, Q < 6 < 7r/2 and the required q is, 



(45) 



(46) 



TT I TT N 
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(47) 
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Figure 5: Probability of success using the qa and q for both algorithms: a. 
< M/N < 1 (left), b. M/N < 1 x IQ-^ (right). 



Fig. I^la shows the probability of success for both algorithms using the 
required number of iterations. We can see from the plot that the minimum 
probability that Grover's algorithm may reach is approximately 50.0% when 
M/N « 0.5 while for the proposed algorithm, the minimum probability is 
87.88% when M/N « 0.2928. Grover's algorithm wiU behave similar to the 
single guess technique for M/N > 0.5 where qc = in that range so that 
Ps = M/N. Although the proposed algorithm is slower than Grover's al- 
gorithm for small M/N by v^, Fig. Elb shows the probability of success for 
both algorithms for small M/N (hard cases where M/N < 1 x 10"^), where we 
can see that the proposed algorithm is more reliable (higher probability) than 
Grover's algorithm. For M/N > 1/3, q = I where the proposed algorithm runs 
in 0(1) to get probability at least 90%, i.e. the problem is easier for multiple 
matches. 



6 Unknown Number of Matches 

In case we do not know the number of matches M in advance, we can apply 
the algorithm shown in ^ for 1 < M < TV by replacing Grover's step with the 
proposed algorithm. The algorithm can be summarised as follows: 

1- Start with m ~ 1 and A ~ 8/7. (where A can take any value between 1 
and 4/3) 

2- Pick an integer j between and ttt, — 1 in a uniform random manner. 

3- Run j iterations of the proposed algorithm on the state: 
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N-1 
1=0 



4- Measure the register and assume i is the output. 

5- If /(i) = 1, then we found a sohition and exit. 

6- Let m = min I Xm, vTV) and go to step 2. 



For the sake of simphcity and to be able to compare the performance of this 
algorithm with that shown in ^ , we will try to follow the same style of analysis 
used in 0]. Before we construct the analysis, we need the following lemmas. 
The first lemma is straightforward using mathematical induction. 

Lemma 6.1. For any positive integer m and real number 9 such that < 6 < 
■k/2, 

771 — 1 

sin^ {{q +l)0)+ sin^ (qe) = m 

q=0 



E- 9 // 1N/,N • 9 / /,^ cos (6') sin (2m6 
sm2 {{q +1)9)+ sm^ (qO) = m ^^iiniO) 



Lemma 6.2. Assume M is the unknown number of matches such that 1 < 
M < N. Let 9 be a real number such that cos (9) ~ I — M/N and < 9 < 7r/2. 
Let m be any positive integer. Let q be any integer picked in a uniform random 
manner between and tti — 1 . Measuring the register after applying q iterations 
of the proposed algorithm starting from the initial state, the probability Pm of 
finding a solution is as follows, 

Pm = ^—^ ( 1 - 



1 / cos (6*) sin {2m9) \ 

1 + cos (9) \ 2TOsin(6l) ) ' 



where, P„ > 0.2725 for m > 1/ sin (9) and small M/N. 

Proof. The average probability of success when applying q iterations of the 
proposed algorithm when < q < m is picked in a uniform random manner is 
as follows, 

m — 1 , , m — 1 

Pn.= E ^P^'^ = ^(i+L(.)) E sin2((g+l)0)+sin2M) 

q=0 q=0 

-1 cos(^) sin(27n^) 



^ l+cos(0) \^ 2msin(e) 

If TO > 1/ sin (9) and A/ < iV then cos {9) w 1, so, 

1 sin (2in9) 1 sin (2ni9) 
''->-2-4^)^-2 T^>0-2725, 

where sin(2m6') < 0.91 for < 6* < tt/2. 
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□ 



We calculate the total expected number of iterations as done in Theorem 3 
in U. Assume that niq > l/sm{9), and Vq = [log;^ m^] . Notice that, niq = 

O (^/nJm\ for 1<M <N, then: 

1- The total expected number of iterations to reach the critical stage, i.e. 
when m > rUq: 






2^-^ 2(A-1) 



rria — i.brjio. 



2- The total expected number of iterations after reaching the critical stage: 

1 °° 1 

- g (0.7275)" A-- = 2(1-0.7275A) -- ^ ^'^-^^ 

The total expected number of iterations whether we reach to the critical 
stage or not is QArUq which is in 0{y/N/M) ioi 1 < M < N. 

When this algorithm employed Grover's algorithm 01, and based on the con- 
dition niG > 1/ sin (26*0) = O (^/nJm\ for 1 < i\/ < 3iV/4, where ma wiU act 
as a lower bound for qo in that range. The total expected number of iterations 
is approximately Smc- For M > 37V/4, niG will increase exponentially where it 
will not be able to approximate qc- Employing the proposed algorithm instead, 
and based on the condition niq > 1/ sin (6) ^ O ( \JN/M J for 1 < A/ < N , the 

total expected number of iterations is approximately 6Amq, i.e. the algorithm 
will be able to handle the whole range, since niq will be able to act as a lower 
bound for q over 1 < M < A^. Fig. El compares between the total expected 
number of iterations for both algorithms taking A = 8/7. 

7 Conclusion 

In Grover's algorithm, the search space is split into two subspaces (the solution 
and non-solution subspaces) then amplifies the amplitudes of the solution states 
by iterating the diffusion operator and the oracle ^2 to find a match with 
high probability in 0{^N/M) for small M/N and in the neighbourhood of 
M/N ~ 1/4 0]. The main idea of using partial diffusion in quantum search 
is to split the subspacc of the solutions into two smaller subspaces. In each 
iteration, one of the solution subspaces will be inverted about the mean (together 
with the non-solution subspace) while the other half will have the sign of their 
amplitudes changed to the negative sign, preparing it to be inverted about the 
mean (together again with the non-solution subspace) in the next iteration. 
The benefit of this alternating inversion is to preserve half the number of the 
solution states at each iteration so as to resist the de- amplification behaviour 
of the standard diffusion operator when reaching the so-called turning points 
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Figure 6: The actual behaviour of the functions representing the total expected 
number of iterations for Grover's algorithm Smc and the proposed algorithm 
GAiriq taking A = 8/7, where the number of iterations is the flooring of the 
values (step function). 



and get the solution with high probability in 0{^ N/M) for \ < M < N. 
Apply the oracle Uf each iteration will switch the entanglement of the two 
solution subspaces with the extra qubit workspace to decide which subspace to 
be inverted about the mean with the non-solution subspace. 

An algorithm for unknown number of matches replacing Grover's step in 
the algorithm shown in 0] is presented, where we showed that the algorithm 
will be able to handle the range 1 < M < iV in 0{\J N/M) compared with 
1 < M < 3iV/4 when using Grover's algorithm. 

We showed that the algorithm will be able to handle the whole possible range 

1 < M < N more reliably using fixed operators in O I yj N/M\ for both known 

and unknown number of matches which makes it more suitable for practical 
purposes. 
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